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Introduction. The launch vehicle (LV) in flight and the dynamic components of loads from the impact of a trapezoidal 
wind gust are considered. It is proposed to determine the dynamic components of the force factors using analytical 
solutions for the structure points accelerations. The work objective is to create a technique for selecting the duration of 
the standard gust, under the influence of which maximum loads are provided in the sections of the LV structure. 
Materials and Methods. The launch vehicle is presented as an uneven beam. The description of its vibrations 1s reduced 
to a system of independent ordinary differential equations that determine the motion of an equivalent system of 
oscillators. The equation of oscillator vibrations under the action of a trapezoidal pulse load is solved by the overlay 
method, and it is reduced to the calculation of the Duhamel integral. It 1s proposed to get the parameters of an 
equivalent system of oscillators based on the results of the calculation of dynamic characteristics for a finite element LV 
model in the Nastran program. 

Results. Analytical relations for the LV structure point accelerations under the action of a trapezoidal wind gust are 
given. For the beam model, test calculations of accelerations were carried out according to the technique proposed in 
this paper. These data are compared to the results of finite element modeling. With the help of analytical solutions, 
dependences are constructed that determine the nature of the change in the magnitude of the bending moment for 
different sections of the launch vehicle when the duration of the wind gust varies. 

Discussion and Conclusions. The presented technique provides building an equivalent dynamic model of systems with 
a large number of degrees of freedom on the example of a LV and obtaining analytical solutions for accelerations of 
points of a mechanical system under trapezoidal external action. These solutions are applicable for the study of dynamic 
loads. The analysis results enable to select the duration of the wind gust, at which maximum loads are reached in the 
sections of the LV structure. Calculations based on the analytical solutions are very economical in terms of time spent. 
They can be used in design calculations for preliminary assessment of loading. 
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Introduction. During the operation of the launch vehicle (LV), loads occur in the elements of its design. We 
are talking about longitudinal and shear forces, bending and torsion moments. Data on these force factors are used for 
strength analysis in the design of new products, experimental development of the design [1], and adaptation of launch 
vehicles for a specific start-up [2]. Loads are divided into quasi-static and dynamic. The quasi-static ones arising in 
flight are calculated from the condition of dynamic equilibrium of the LV as a solid body, taking into account the 
permissible parameters of the ascent trajectory. 

Generalized beam models are usually used to calculate dynamic loads. Such loads are determined by the results 
of solving the equation of the elastic LV motion, which in general is a partial differential equation. Methods based on 
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the application of decomposition of the solution by the eigenvibration tones of the structure are able to provide a high 
speed of calculation in combination with sufficient accuracy of the results [3]. As shown in [4], using the method of 
decomposition by the forms of natural oscillations, it is possible to proceed to a system of independent ordinary 
differential equations. They describe: 

— centroidal motion of the LV, 

— rotation of the longitudinal axis of the LV relative to the center of mass, 

— transverse elastic modes of the LV in flight. 

The transition to independent equations describing elastic modes of the LV means that the distributed 
parameters of oscillators will be taken into account, each of which is a single-degree-of-freedom system. The motion of 
such an oscillator can be considered independently of the others, and a solution can be obtained for each of them using 
well-known methods of oscillation theory. 

To calculate the loads in flight, not the entire ascent trajectory is considered, but only some of its points, the so- 
called load cases, characterized by the extreme value of individual parameters affecting loading, or the maximum value 
of loads on individual structural elements. One of the most important cases of loading 1s the atmospheric flight of the 
LV [1, 5]. The influence of a turbulent atmosphere on the LV loading can be determined by statistical methods [6, 7], or 
within the framework of a conservative approach, when the maximum possible (with some level of probability) wind 
characteristics are taken into account. This paper discusses the second approach. A single specified wind gust is 
accepted as an external dynamic effect. The profile of the specified gust, which characterizes the change in wind speed 
over time, can be set in a trapezoidal [8], cosine, or sinusoidal form [9]. In this paper, we will consider the launch 
vehicle motion under the trapezoidal gust action. The duration of the specified gust is usually selected to be comparable 
with the period of the lowest transverse tone of the LV’ vibrations. At the same time, there are often demands for its 
variation to achieve maximum efforts in the LV sections [8, 10]. The complexity of calculations using standard finite 
element (FE) analysis programs is due to: 

— the need to vary the parameters of external action, 

— a large number of calculated cases, 

—a variety of design options and configurations of structures at the stage of design calculations [9]. 

This paper objective is to develop a methodology for selecting the duration of the specified gust using 
analytical solutions obtained for a simplified dynamic LV model presented as an equivalent system of oscillators. 

A semi-analytical approach using the Duhamel integral was successfully applied in [11] for hydroelastic 
analysis of ships. In [12] and a number of other works, the Duhamel integral is used as part of the problem solution of 
loading bridges with moving loads. In this paper, the Duhamel integral is used for analytical solutions to the LV 
reaction to the short-term impact of a wind gust in flight. 

Materials and Methods. At the stages of preliminary design, it is advisable to use flat design schemes for 
beam models. With the simplicity and speed of the solution, they enable to determine the motion variables and internal 
forces (with an accuracy acceptable for this stage of design) [13]. Imagine the LV in the form of an elastic beam with 
variable length mass and stiffness. Assume the usual assumptions for the resistance of materials, including the 
hypothesis for the smallness of elastic deformations. To determine the internal force factors in the LV section, we use 
the acceleration (overload) method [1, 4], which can be interpreted as a cross-section method adapted for dynamic 
calculation. In this case, internal forces are found from the conditions of static equilibrium of mentally cut off parts of 
the structure under the action of external distributed loads, supplemented by D'Alembert’s forces of inertia, and the 
desired internal forces. Quasi-static and dynamic values of force factors are determined separately based on pre- 
calculated accelerations, and then summed up [4]. 

This paper considers the issue of determining the dynamic loading of the LV in the transverse direction under 
the action of a wind gust, whose speed is directed perpendicular to the longitudinal axis of the LV. It is assumed that the 
loading in the longitudinal direction can be calculated independently. It is not considered in this paper. 

To determine dynamic accelerations, the LV is presented as a free elastic beam. Its motion is studied in the 
vicinity of the moment of time corresponding to the load case under consideration, and is described in deviations from 
the state of dynamic equilibrium in which the LV was before the wind gust, moving along the nominal (undisturbed) 
trajectory. Here, such parameters as the mass and moment of inertia of the LV, the motor power, and the angle of 
projection are assumed to be constant and equal to the characteristics of the considered point of the nominal trajectory. 
The perturbed motion of the elastic LV is investigated in a fixed coordinate system associated with the position that the 
LV occupied at the time of the calculation. The perturbed motion will be a combination of plane-parallel motion of the 
LV as a rigid body in the plane in which the dynamic load is applied, and elastic modes of the structure. We do not take 
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into account the control system response, i.e., for the stabilization machine, we assume a long delay time compared to 
the time of application of the dynamic load. In general, under the influence of a wind gust, together with elastic modes 
of the body, the LV starts to move as a rigid body. The transverse component of the aerodynamic force, which is 
considered proportional to it at a small angle of attack, changes its value due to the displacement of the structure in the 
direction of the wind gust and rotation relative to the forward flow vector. The projection of gravity on the transverse 
axis of the LV also changes. To properly account for these changes, the equations of LV motion must be integrated with 
the equations describing the logic of the automatic stabilization, which is impossible at the early stages of design. 
Taking into account the significant mass and the moment of inertia of the LV, we will consider small: 

— the angle of rotation of the LV as a rigid body for the time of calculation; 

— the rate of displacement of the LV center of mass in the direction of the wind. 

This enables to ignore the impact of the above changes. The angle-of-attack increment (and, consequently, the 
transverse component of the aerodynamic force) is considered to depend only on the magnitude of the wind gust speed 
given as a function of time. Thus, taking into account the accepted assumptions, the aerodynamic load in the transverse 
direction is a load distributed along the length of the beam with a time-dependent proportionality coefficient. The law of 
distribution of aerodynamic load along the length of the LV is determined experimentally and is considered to be 
known in advance. The law of change of the proportionality coefficient (angle of attack) from time to time is 
determined by the selection of the profile of the specified wind gust. 

The motion of the LV modeled as an elastic beam can be described using the well-known equation of forced 
transverse vibrations of the beam, written with the account 2 Voigt's hypothesis: 


m(x) 229 + (14+ h2)F [B@) 2] = a0), (1) 


where m(x) — mass per unit length; B(x) — bending stiffness; q(x,t) — distributed external load; h — friction 








factor. 
This equation should be supplemented with a boundary condition: the internal forces in the initial and final 
sections are zero. It means: 


2 (BG) 2) = 0, B(x) 2E9 = 0, mupyx=0,x =L, (2) 
where L — the LV length. 


In this paper, the aerodynamic force is taken as an external distributed load, which can be presented as a 
product of functions: 


q(x,t) = R(t)Y,(%), (3) 
where R(t) — the function that determines the temporal variability of the aerodynamic force and varies according to 
the trapezoidal law in accordance with the wind gust model adopted in this paper; Y,(x) — the function of the 


aerodynamic force distribution along the LV length. 

Let us consider free LV vibrations without taking into account friction forces (at q(x,t) = 0, h=0). We 
substitute the variable separation method y(x,t) = f(x):q(t). In this case, from equation (1) with boundary 
conditions (2), it is possible to arrive at an Rae differential equation 


*[B@) 2] - p?ma@f@) = 0, (4) 





with boundary conditions: 








S[B@) LO] = 0, Bo) LO =o, npux=0,x=1 (5) 


Solution (4) with conditions (5) is a classical Sturm-Liouville problem. Solving it, one can find a set of 
eigenforms f;(x) and eigenfrequencies p; of the beam in question Vj = 1,2,...). It is known’ that some solutions to 
system (4) correspond to zero natural frequencies. The forms corresponding to zero natural frequencies determine the 
translational motion of the LV as a rigid body together with the center of mass and rotation around the center of mass: 
f-1 = 1, fo = x — Xc, where x; — coordinate of the LV center of mass. 

It should be noted that the LV mass and stiffness characteristics most often have a piecewise constant pattern 
of distribution. In this case, the equations of form (1) and (4) should be written separately for each homogeneous beam 
section with boundary conditions at the junctions of the sections, as in the derivation of the ratios of the initial 
parameters method [1,14]. In the given paper, this entry is omitted, because the calculation of the dynamic 
characteristics (modal analysis) of structures is carried out numerically, using the finite element method. 

Let us imagine the forced vibrations of an elastic beam modeling the LV structure in the form of 
decomposition according to the natural modes. Assume that the beam stiffness axis passes through its center of mass. 
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To move the points of the LV axis, we write: 

y(x,t) = ye(t) + 9) — xc) + D8 HR) “aj ©, (6) 
where y(t) — displacement of the beam center of mass; 9(t) — angle of rotation of the axis of the undeformed beam; 
f; — beam eigenmode corresponding to tone numbered j; q,;(t) — generalized coordinate corresponding to tone 
numbered j; N — number of elastic tones taken into account. 

After substituting (6) into (1) and applying the Bubnov-Galerkin procedure, we can arrive at an ordinary 
differential system with constant coefficients: 
MYc(t) = Qy, 
19(t) = Qo. (7) 
m,(q;(t) + 2njq,(t) + p74q;(t)) = Q(t) G = 1,2,...,N) 
Here, m — LV mass; J — moment of inertia relative to the axis passing through the LV center of mass perpendicular to 
the plane of rotation; m; — reduced (generalized) mass for the j-th vibration tone and determined from the formula: 


L 
mj = J, m(x)f/ (x) dx. 
The generalized forces in expression (7), taking into account (3), are defined as follows: 


i 
Q, = R(t) i Y, (x) dx = RU); (8) 
L 

Qs = R(t) [. Yq (x) (« — xX¢)dx = R(t)Mg, (9) 

i 
Q; = R(t) fe Y,(x) fjdx = R(t)Qo;. (10) 
Here, Q, — maximum value of the main transverse aerodynamic load vector; M, — maximum value of the main 
transverse aerodynamic moment reduced to the LV center of mass; Q); —- maximum value of the generalized force 


corresponding to the generalized coordinate qj. 

The first two equations in (7) define the law of change of accelerations of the LV points in the process of 
translational and rotational motion of the LV as a rigid body. The last equation in (7) defines the law of motion of an 
equivalent system of oscillators. 

Consider the motion of one oscillator under the action of a trapezoidal external load, which: 

— increases from zero to Qp over time 6, 

— maintains a constant value over time 0, 

— drops to zero over time 6. 

For convenience, we omit the indices characterizing the tone number. Then we will rewrite the differential 
equation of the oscillator motion taking into account (10) in the form: 


G+ 2ng + p?q = =“ RO), (11) 
where m — oscillator mass; p — angular frequency of natural oscillations, expressed in radians per second; n — 
damping coefficient (determines the oscillator damping); Q(t) = Q)R(t) — law of variation of external load. 

We represent function R(t) as a set of four linear functions: 

Re) =o. 2 R(t) = —8"'(t — 9), 
R3(t) = —6°-*(t — 6-8), R,(t) = 5 *(t — 26 — 9). 

Accordingly, external load Q(t) is a combination of four linear loads Q;(t) = Q)R;(t), (i = 1, 2,3,4). Load 
Q,(t) is applied from moment t = 0; Q,(t) — from moment t = 6; Q3(t) — from moment t = 6 + 6; Q,(t) — from 
moment t = 8 + 26. 

To determine the system response to external actions, we divide the entire duration of the load into four 
intervals (Fig. 1). 


Q(t),N 





' Ist section | 2nd section : 3rd section; 4th section 


Fig. 1. Diagram of the external trapezoidal load 
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In accordance with the superimposition method’, let us represent the response of the linear system under 
consideration to an external action in the form of the Duhamel’s integral as the sum of responses to a set of 


independently applied elementary impulses: 
t 


a(t) = | P(e- 8) (0)ab, 
0 
Here, P(t — 8) — the shifted time law of variation of external action, and Y(9) characterizes the system 


response to a unit impulse input. The response of mechanical system q,(t) to a linearly increasing load P(t) = kt can 


be expressed through the system response to unit impulse Y(9) and to a suddenly applied single load Y, (9): 
t 


q(t) =k | Y,(9)d8, Y,(t) = | Y(9)d9. 
0 


0 


The response of a one-degree-of-freedom mechanical system with damping coefficient n to a unit pulse will 


have the form’: 


1 
Y(9) = —————-e™”* sin ( p2 — n29), 


We introduce a notation for the frequency of damped vibrations p, = ,/p* — n2. Let us calculate the system 


response to a suddenly applied unit load: 
t 








1 n 
Viiv) = | Y(9)d9 = E —e (cos( t) + —sin( ‘))} 
1 mp? P41 D, P41 
The response to a linearly increasing load: 
_ t _  k (2n =i at (an? —nt o: 
q(t) =k in Y,(9)d9 = aa (s [—1 +e cos(p,t)] Les ( 2 1) e~™ sin(p,t) + th (12) 


We accept the notations: A,;= t — 6, A,= A, — 98, A3= A, — 4. 


We introduce a function containing the harmonic terms of solution (12): 
n he re yh 
d(9) = 2= (cos(p, -0) +—sin(p, -9)) e ™ ——sin(p,:8)e-™. 
p P41 P4 


We denote by q,, the motion under the action of a force statically applied to the system Q) = g,,mp” and take 
into account that k = Q,)/6. The total movement of the oscillator under the action of a combination of loads Q;(t) at 
every time point will be the sum of the corresponding solutions (12). The system response described by equation (11) to 


the external impact of the trapezoidal profile will have the form: 


daS "ld = 2n/p> +t] at0<t <6, 
oe derS *[d(t) — d(A,) + §] at8<t<6+6, 13) 
ded 1 [d(t) — d(A,) — d(A,) + 2n/p?—A,] atO9+5<t< 0+ 26, 
der *(d(t) — d(A,) — d(A2) + d(A3)) att > 0 + 26. 


The law of variation of the oscillator accelerations can be obtained from double differentiation in time of 


d2 


aa (d (9)). We differentiate and record the result taking into account 


expression (13). Let us introduce function g(9) = 


the number of the oscillation tone: 


Nn; Nn; 1 
g;(8) = G = pt ;) (24 (costo, -9)+ —-sin(p, ; )) = —sin(p, ; )) — 
Dj P41 j P41 j 
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=on; (22 (py, sin(p, ; -9) ‘Pad; cos(p, ; -9)) = cos(p, ; -8))| ew. (14) 


We rewrite the expression for the frequency of damped natural vibrations: 


Pij= Ip; ie (15) 


To determine the accelerations of the LV points, we differentiate twice (6): 


W(t, x) = ¥e(t) + OO) — xe) + VL HOG (O. (16) 
Let us express the accelerations of the generalized coordinates from the first two equations of system (7) taking 


into account (8) and (9). We represent the accelerations of the points of the LV axis that it acquires through moving as a 


rigid body: 
a(t) = ¥c(t) + 8(E) (x — xc) = R(t)(Qa/m+t M,/1). (17) 
In addition, we take into account the law of variation of function R(t) from time: 
5 1t atO <t <6, 
R(t) = 1 até6<t<68+46, 
—§ 4(t — 26-6) at0+6<t<0+4+26, 
0 att >0+ 26. 
We get: 
5-'t(Q,/m+ M,/1) ac 0=t.< 5, 
ioe Q,/m+M,/1 atOé<t<0+46, (18) 
—5-1(t — 26-6)(Q,/m+M,/I) w8+6<5t<0+ 26, 
0 att. 0+ 20. 


For accelerations of generalized coordinates corresponding to the elastic vibration tones, as a result of double 
differentiation of expression (13), taking into account (14), we obtain: 
5 der 59; (t) at0<t <6, 
8 qer j (gj(t) — 9;(Ay)) até6<t<6+6, 
8 der j(gj(t) — gj(Ar) — gj(A2)) at +5 <b <0 +26, 
8 der j (gj(t) — 9j(Ar) — gj(A2) + gj(A3)) att > 0+ 28 


Taking into account (16) and (17), the law of variation in the accelerations of the elastic LV axis points under 


qj(t) = (19) 


the action of a trapezoidal wind gust will have the form: 
VE, x) = a(t) + Yj ffO)GO, (20) 
where a(t) is determined from expression (18), g;(t) — from expression (19). 
While investigating the LV elastic modes, damping is traditionally taken on the basis of the data obtained from 
the results of full-scale dynamic tests and presented in the form of values of logarithmic decrements D;. Then, the 
coefficient determining the damping parameter and included in expressions (14) and (15) can be calculated from 


formula: 
The change in generalized coordinate q; included in formula (19) under the action of a statically applied 
generalized force Qo; is determined from expression: 


4 Qoj 4 fy Yai) ax 
Pp Mj = PF Jo mx) F? (x) dx 


cr J ~ (22) 


Taking into account (14), (15), (18), (19), (21), (22), formula (20) is an analytical expression that defines the 
functions of the acceleration change over time for the points of the LV axis under the action of an external transverse 
aerodynamic force varying according to the trapezoidal law, taking into account the impact of dissipative forces. 


Knowing the law of acceleration variation, it 1s possible to determine the dynamic, and then the total structure 
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loads acting in the LV sections using the known methods [4]. We calculate the bending moment due to inertial forces 
from the elastic modes of the LV structure: 
M(x,t) = — Vj Mj (1) GO. (23) 
Here, M;(x) — function of distributing a unit (with acceleration q;(t), equal to one) bending moment for the j-th 
vibration tone along the LV length. It can be found from formula: 
M,(x) = fr fo m(x)f; (x)dx dx. (24) 

As noted above, for the modal analysis in this work, the FE method was used. This approach is due to the fact 
that in practice, the dynamic LV model has rather complicated structure. It includes substructures, and their own 
dynamics cannot be neglected. Substructures can be attached to the LV body in one section, or be located parallel to the 
longitudinal axis of the LV and have several attachment points. In this case, the calculation of dynamic characteristics in 
a continuum setting is a complex mathematical problem. In addition, dynamic models of individual substructures are 
presented by development companies in a condensed (matrix) form in the Nastran format. For this reason, it will be 
optimal to use the Nastran engineering analysis software package to calculate the dynamic characteristics of the 
structure. 

However, the use of standard FE analysis programs for the calculation of dynamic loading is fraught with 
certain difficulties. These include the need to pre-construct an equivalent model of the external aerodynamic load 
suitable for use in the FE analysis program [15], which in itself is quite difficult. In addition, there are difficulties 
associated with processing the calculation results. The use of the postprocessor functionality for analyzing the results is 
extremely time-consuming and requires a large number of manual operations. Another way involves the application of 
additional software for processing large array of numeric data”. In this paper, an approach is proposed in which the FE 
analysis program is used only for modal analysis. In this case, dynamic loading is calculated using specially developed 
software that allows you to vary external loads and automatically process the calculation results. 

For a complex LV design, equations (7) retain their form [3]. The standard output information of the Nastran 
program can be the basis for obtaining the parameters of an equivalent system of oscillators, LV mass and moment of 
inertia, as well as for calculating the generalized forces included in the third equation of system (7). To form the left 
side of the third equation in (7), the values of natural frequencies p, (Radians) and generalized masses m, (Generalized 
mass) are required. To determine the generalized forces in the right part of the third equation in (7), the eigenmode 
functions f;(x) (Eigenvector) are needed. When calculating dynamic inertial loads according to (23), instead of (24), it 
is more convenient to use unit inertial loads (forces and moments), which are output by the Nastran program after the 
standard application of forces when calculating eigenforms and frequencies. Unit inertial loads are output separately for 
each vibration tone and multiplied by the square of the natural frequency, which should be taken into account for their 
correct use. 

Research Results. To carry out test calculations, a medium class tandem launcher is considered. Figure 2 
shows the type of functions of linear bending stiffness B(x) and mass m(x), as well as the distribution of concentrated 


MASSES Megcp(X) = Meocp A(X — x,) along the length of the considered LV (the Dirac delta function is denoted by A). 
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m(x), kg/m 


Mecocp (x) ) kg 


kg-N 
B(x), =— 


c) 


Fig. 2. Distribution of mass and stiffness characteristics along the LV length: 
linear mass (a), concentrated mass (b), bending stiffness (c) 


To check the obtained analytical solutions under the impact of a trapezoidal external load, the dynamic LV 
accelerations are calculated according to formula (20), taking into account (14), (15), (18), (19), (21), (22). At the same 
time, 5 elastic tones of natural transverse LV vibrations are considered. Characteristics p;, m, and f;(x) were obtained 
from the calculation results in the MSC Nastran software package using the solution sequence for modal analysis of 
natural vibrations (SOL 103). The LV dynamic finite element model is presented as a set of beam elements with various 
inertial and stiffness characteristics. The following elements are elastically or rigidly attached to the beam elements: 

— elements describing the inertial properties of devices, aggregates, parts of the block structure; 

— condensed models of individual blocks presented in the digital matrix form. 

The finite element LV model includes about 1000 components. The elements simulating restraints were not 
used to preserve the ability of the LV to move as a rigid body. 

According to the same finite element LV model in the MSC Nastran software package, the accelerations are 
calculated through the SOL 119 solution sequence used for modal transient analysis. At the same time, all tones of 
natural vibrations in the range up to 100 Hz were taken into account in the modal decomposition. The aerodynamic load 
is represented by transverse linear loads distributed over all beam elements simulating the LV structure. In addition, the 
trapezoidal law of the aerodynamic load variation over time 1s given. 

Figure 3 shows the results of a comparative analysis of accelerations obtained using two different approaches 
A(t) = V(t,xX9) of a certain point of the LV axis with coordinate x = x9. It can be seen that the following two 


solutions agree well: 
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— the solution obtained for accelerations using a simplified LV model based on the analytical relations given 


in this paper; 
— numerical solution obtained from the complete finite element LV model. 
a,(t), 30.0 
2 
m/s" 20.0 


10.0 
0.0 + 
-10.0 
-20.0 
-30.0 + 





Fig. 3. Comparison of transverse accelerations obtained through two different 
— finite element modeling, ® °° °° — analytical solution for a simplified model 





approaches: 

The peak value of the bending moment selected from the time process determines the level of equivalent forces 
taken to carry out the strength calculation, and acts as a variable parameter when changing the parameters of external 
action in the transverse direction [16]. 

Dependences are obtained based on the values of dynamic accelerations (Fig. 4). They show how the change in 
peak bending moment M or different LV sections are related to parameter 8, which characterizes the gust duration. In 
Figure 4, the bending moment values are presented in the form of dimensionless quantities M*. They are calculated 
through dividing the dimensional bending moment by the maximum value for a given section (e.g., for section x = 0.3L 


maximum value M,,,,, = 4.6: 10° N-m), found when value 0 varies over the entire range under consideration. 


M" 1.0 
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Fig. 4. Dependence of the value of dimensionless bending moment M* on the duration of wind gust 0*, expressed in fractions of 
period T, for various LV sections: zone 1 (x = 0.0 — 0.2L) (a); zone 2 (x = 0.2L — 0.4L) (6); 
zone 3 (x = 0.4L — 0.55L) (c); zone 4 (x = 0.55L — 0.75L) (d); 
zone 5 (x = 0.75L — 0.9L) (e); zone 6 (x = 0.9L — L) (f) 


So, on the graphs, the duration of the gust action is represented by dimensionless value 0° obtained through 
dividing parameter @ by period T, of the first tone of the LV vibrations. Each line of the graph corresponds to one 
section of the LV. All LV sections are grouped according to the nature of function M*(6*) and are shown on various 


graphs, and the LV length is appropriately divided into zones (Fig. 5). 






— 


Fig. 5. LV cross-section zones 





It can be seen from Figure 4 that for the first zone, the maximum bending moment is reached already at 0 
values at least 15% of period T, of the first vibration tone. For the second and third zones — slightly more than half of 
period T,. A further increase in 0 values does not affect the magnitudes of the maximum values of the bending moment. 
For zones 4-6, the maximum value of the bending moment turns out to be local and is within the zone of 8 values close 
to the value of half of period T, of the first tone of the LV vibrations. The results obtained fully correspond to the results 


of finite element modeling carried out earlier®. 


° Malykhina OI, Glugovskii MS. Analiz vliyaniya profilya poryva vetra na velichinu korpusnykh nagruzok rakety-nositelya v poletnykh sluchayakh 
nagruzhentya. In: Proc. VI All-Russian Sci.-Tech. Conf. Samara. 2019;1:133-138. (In Russ.) 
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Thus, to obtain the maximum values of the bending moment in the sections of the considered LV, a dynamic 
analysis of the behavior of the structure with external action in the form of a wind gust of is required. Its duration is 
determined by parameter 9, close in value to half of the period of the first tone of the LV vibrations. 

Discussion and Conclusions. Using the superposition method, analytical solutions are obtained that describe 
the motion of a one-degree-of-freedom system, which is affected by the friction force and an external force varying 
according to the trapezoidal law. The method of application of the obtained analytical solutions for systems with many 
degrees of freedom is given. A good coincidence of two types of solutions is shown: 

— analytical one, for accelerations of LV points found from a simplified model; 

— numerical data obtained from the complete finite element LV model. 

It is shown that analytical solutions can be used to analyze dynamic force factors to select the duration of a 
wind gust, under the impact of which maximum loads are achieved in the sections of the LV structure. Similarly, it is 
possible to analyze the overloads that are achieved in the LV sections (e.g., at the installation points of measurement 
systems). 

In addition, the proposed methodology provides building a full cycle of the load analysis pre-calculation in the 
case when an analytical representation of the external dynamic load is possible. The load analysis based on analytical 
solutions 1s very economical in terms of calculation time, and it can be a remarkable alternative to finite element 
modeling at the design stage, when a large number of combinations of external loads and configurations of the design 
under development are studied. Finite element analysis of the detailed model in this case can be used as a refine final 


calculation. 
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